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Abstract 

We study the vacuum statistics of ensembles of M theory compactifications on G 2 
holonomy manifolds with fluxes, and of ensembles of Freund-Rubin vacua. We discuss 
similarities and differences between these and Type IIB flux landscapes. For the G 2 
ensembles, we find that large volumes are strongly suppressed, and for both, unlike the IIB 
case, the distribution of cosmological constants is non-uniform. We also argue that these 
ensembles typically have exponentially more non-supersymmetric than supersymmetric 
vacua, and show that supersymmetry is virtually always broken at a high scale. 
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1 Introduction and Summary 



String theory has long held the promise to provide us with a complete and final description 
of the laws of physics in our universe. Yet from the start it was clear that the theory 
allows many vacua we clearly do not live in — ten dimensional Minkowski space being the 
simplest example. Over the past two decades, the existence of an infinite set of distinct 
stable vacua in which we obviously do not reside was firmly established. For some of 
those, such as Type IIB string theory on AdS 5 x S 5 with flux, we even have complete 
nonperturbative descriptions. The vast majority of these well understood vacua have 
unbroken supersymmetry and come either with massless scalars or a negative cosmological 
constant of the order of the KK scale, all properties we certainly do not observe. Moreover, 
even though constructing vacua with light charged particles is fairly straightforward, 
typically such vacua are very different from the standard model at low energies. 

The implicit assumption, of course, was that at least one vacuum should exist with 
broken supersymmetry, with just the standard model at low energies, and a couple of 
rather astonishing hierarchies between the cosmological, electroweak and Planck scales. 
One could hope that some deep dynamical mechanism uniquely selects such a vacuum, 
with all its peculiarities, rather than one from the enormous set of consistent vacua already 
firmly established, or from the presumably even much more enormous set of vacua that do 
not satisfy the rather artificial constraint of being under our full computational control at 
the present stage of theoretical development. Thus far however, no indication whatsoever 
has emerged that this is indeed the case, and anyway there is no realistic hope that the 
"correct" vacuum could be computed from such hypothetical dynamical principles with 
what we currently know about string theory. 

It thus seems more modest and productive to follow the historically successful approach 
of trying to construct models that are compatible both with experimental and theoretical 
constraints. Whilst progress is continually being made in the program of constructing 
string theory vacua which meet rough observational requirements such as obtaining the 
standard model spectrum at low energies, no massless scalars or a small positive cosmo- 
logical constant, no explicit model has been constructed combining all of these features. 
Clearly though, the examples studied thus far represent only an infinitesimal part of the 
tip of the iceberg, as computational complexity quickly turns into a major obstacle. The 
iceberg in question has been called the string theory landscape, and the problem we face 
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is to find places in this landscape that could describe our universe. 

Before embarking on this program however, it would evidently be useful to know what 
our chances of success are. For example it could save us a lot of time and effort if we man- 
aged to exclude large regions of the landscape, or more precisely large classes of models, 
before attempting to find realistic string theory vacua there. Traditionally the approach 
has been to focus on the most easily controllable constraints such as light charged particle 
spectra, while ignoring issues such as moduli stabilization, supersymmetry breaking or 
the cosmological constant. However these constraints are equally important, and it would 
be very useful to estimate how much they reduce the set of possibilities. For instance, 
one could propose a certain compactification with just the standard model living on a set 
of branes, and with very large compactification volume to explain the hierarchy between 
the electroweak and Planck scales, along the lines of PQ|. This would already seem like a 
considerable triumph. A courageous model builder might then go on and try to turn on 
various fluxes to stabilize all moduli within the required region of moduli space, following 
the scenario of [2] or [3]. But this theorist will plausibly end his project disillusioned, 
finding not a single vacuum with volume even remotely close to what is required for this 
model to work. Indeed, it is generally true in these scenarios that flux vacua rapidly be- 
come more scarce at larger compactification volumes, and completely cease to exist even 
at very modest sizes. Clearly, our model builder would have benefited greatly from some 
simple estimates of distributions of actual vacua over the parameter space of the model. 

In it was furthermore pointed out that a sufficiently fine "discretuum" of string 
theory vacua could accommodate an extremely small but nonzero cosmological constant. 
To find out whether this idea can be realized in a given ensemble of vacua, and more 
generally to analyze how strongly constraints on various parameters reduce the number 
of possibilities within a given class of models (possibly down to zero as in the example 
above), one would like to have an estimate of how many vacua with certain properties lie in 
a given region of the landscape. In other words, one needs to study the statistics of vacua 
in this region (where "statistics" does not refer to any probability measure, but simply 
to number distributions on parameter space). This program was initiated in j3] and has 
been developed, thus far mainly for Type IIB flux vacua, in |H1 IZ1 El 13 HH1 HH H21 HSl HH 
A number of ideas on the important question of distributions of supersymmetry 
breaking scales were proposed in [T7| [TBI HH] , and subsequently significantly qualified in 
|20j . Further considerations on this topic have been made from an effective field theory 
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point of view in [2UI22]- Distributions of gauge groups and particle spectra were studied 
in detail in [23 121] • A number of effective field theory models of the landscape and 
phenomenological consequences have been analyzed in [221 I2IH 123 12H] , and some aspects 
of the statistics of string vacua beyond critical theories have been explored in [29]. Finally, 
more critical viewpoints of the string landscape have appeared in pHH ETT] . 

One class of vacua that are interesting from a phenomenological point of view are M 
theory compactifications on manifolds of Gi holonomy 32] . Because of string duality 
there are presumably at least as many such vacua as there are heterotic string vacua on 
Calabi-Yau threefolds. It is thus important to get an idea of how moduli, compactifica- 
tion volumes, cosmological constants and supersymmetry breaking scales are distributed 
within this class of models. This will be the subject of the main part of this paper. 

Although guidance for model building already provides enough motivation for statisti- 
cal studies, one could be more ambitious and try to infer general properties of distributions 
of string theory vacua from statistical analysis of particular regions of the landscape com- 
bined with genericity arguments. To get some confidence in such results, one would need 
to explore and compare as many regions of the landscape as possible. Thus far how- 
ever, only Type IIB flux vacua have been analyzed in detail. Because of string duality, 
one might optimistically hope that such IIB flux vacua could be representative — in the 
sense that they constitute a significant fraction of all string vacua, perhaps constrained to 
have some additional properties such as supersymmetry in the UV. However, it should be 
noted that the studies in [S| - JE] are limited to vacua described as Calabi-Yau orientifolds 
at moderately large volume and moderately weak string coupling, and strictly speaking 
these represent only a corner of the string theory landscape. In principle, distributions of 
observables could change dramatically as one explores different regions of the landscape. 
It was suggested in [33] that the set of four dimensional string and M theory vacua with 
M = 1 or no supersymmetry is a disconnected space whose different components repre- 
sent qualitatively different low energy physics, and this could translate into very different 
statistical properties. This gives us additional motivation to study the statistics of G 2 
vacua and compare to the IIB case. 

Another branch of the landscape whose statistical analysis we initiate here is the 
set of Freund- Rubin vacua [31], i.e. M theory compactifications on Einstein manifolds 
with positive scalar curvature. Their properties are quite different from more familiar 
compactifications on special holonomy manifolds, and as we will see, this is reflected in 
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their vacuum distributions which are very different as well. 

Finally, one could take things one step further and altogether discard the idea that 
some mysterious dynamical mechanism has uniquely selected our vacuum and in particu- 
lar picked the bizarre and exceedingly nongeneric scale hierarchies which just happen to 
be also necessary to make structure formation and atoms other than hydrogen or helium 
possible (S3 E3 EE EE] • Instead, one could start with the hypothesis that a "multi- 
verse" exists in which a huge number of vacua is actually realized, and in particular that 
we observe ourselves to be in a vacuum with such large scale hierarchies simply because 
this is needed for structure and atoms, and therefore observers, to exist. To make direct 
predictions from just string theory in such a framework, one would need to know the 
probability measure on at least the part of the landscape compatible with a number of 
basic requirements. There is no established way of computing these probabilities at this 
time, but as an additional working hypothesis one might assign for example in a given 
ensemble of flux vacua approximately the same probability to every choice of flux. One 
could refine this by restricting this equal probability postulate to subsets of vacua with 
fixed values of parameters relevant for cosmology, such as the vacuum energy, as one does 
for microstates with equal energy in the microcanonical ensemble of statistical mechanics. 
Different choices of these relevant parameters might then be weighted by cosmological 
considerations (up to the extent that this is needed, as some will be effectively fixed by 
environmental requirements). Under such hypotheses, suitable number distributions can 
be interpreted as probability distributions, and one can test the hypotheses that went in 
by Bayesian inference. 

It should be emphasized that this framework is significantly more predictive than the 
traditional model building approach of simply considering any model compatible with 
current observations. Under these simple hypotheses, the landscape picture together with 
a few rough environmental principles gives a new notion of naturalness for effective field 
theories, which translates into a set of rules for model building. This turns out to lead 
to very distinct models which do not need contrived engineering to fit known data, and 
which moreover give very specific predictions, including many unambiguous signatures at 

lhc mum 123. 

More theoretical data on number distributions obtained from string theory would 
obviously be very useful to make further progress in this area, and although as we discussed 
there are several other motivations for studying the statistics of string and M theory vacua, 
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we consider this to be a very important one. 

The analysis of the M theory vacua we study in this paper is somewhat technical, so 
for the ease of the reader, in the remainder of this introduction, after a brief review of the 
Type IIB flux landscape, we present a summary of our main results. 

1.1 Review of Type IIB statistics 

For comparison, let us briefly review what is known about the statistics of IIB flux vacua 
of Calabi-Yau orientifolds. There is a natural splitting of Kahler and complex structure 
moduli in this context. Turning on flux induces a superpotential jlU] which only depends 
on the complex structure moduli. In suitable circumstances the Kahler moduli can be 
stabilized by nonperturbative effects, along the lines of It is reasonable to ignore 

the Kahler moduli altogether as far as vacuum statistics is concerned, because (1) the main 
contribution to vacuum multiplicities comes from the huge number of different fluxes, (2) 
at sufficiently large volume the Kahler moduli do not influence the positions in complex 
structure moduli space significantly, and (3) practically, the Kahler sector of is less under 
control and more difficult to treat systematically. 

The results for the distributions of moduli, cosmlogical constants, volumes and super- 
symmetry breaking scales for IIB flux vacua are roughly as follows. 

1. The number of supersymmetric vacua in a region 1Z of complex structure moduli 
space is estimated by 

K, c in) = (2v ff )fe3 / detOR + wl) (1.1) 

where L is the induced D3 brane charge by 03, 07 and D7 branes, b% is the third 
Betti number of the Calabi-Yau manifold, R is the curvature form on the moduli 
space and u the Kahler form. Actually this expression gives an index rather than 
an absolute number — it counts vacua with signs, so it is strictly speaking a lower 
bound. Essentially this expression implies that vacua are uniformly distributed 
over moduli space, except when the curvature part becomes important, which is 
the case near conifold degenerations. The above result was verified by Monte Carlo 
experiments in [HH FTf] . 
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2. The cosmological constant for supersymmetric vacua is A = — |W| 2 . 4 Its distri- 
bution for values much smaller than the string scale was found in 9 J to be essentially 
uniform, i.e. 

dj\f ~ Atot dA. (1.2) 
Here Atot ~ L b:i /b$\ is the total number of flux vacua. 

3. The compactification volume V is stabilized by nonperturbative D3-instanton effects 
and/or gaugino condensates, both of which give contributions ~ e ~cv 2 ^ ^ Q ^ e 
superpotential, where c < 1 decreases when 63 increases. In the scenario of [2|, these 
have to balance against the contribution Wq from the fluxes. That is, at sufficiently 
large V (or equivalently sufficiently small Wq), e ~ cV2/3 ^ \Y . Since | W 1 2 ~ A is 
uniformly distributed according to p. 2)1 . this gives the volume distribution 

dN~N tot e- 2cV2IZ d{V 2 /% (1.3) 

Large volumes are therefore exponentially suppressed, and for reasonable values of 
L and 63, the maximal volume will be of the order V max ~ ( log ^" tot ) 3 ^ 2 ~ (6 3 /c) 3 / 2 . 5 
One could also interpret the IIB complex structures fixed by the fluxes to be mirror 
to IIA Kahler moduli. Then the distribution of IIA compactification volumes can 
be shown to be jH] 

dN~^d{V n T) (1-4) 

for Vjja S> 1. Here k is a constant weakly decreasing with increasing 63. Again, 
large volumes are suppressed, now bounded by Vjja < (ekL/b 3 ) 6 . 

4. The flux potential has nonsupersymmetric minima as well. The supersymmetry 
breaking scale F = e K/2 \DW\ = M s 2 usy for F < 1 is distributed as [20] 

dAf~Af tot dF (1.5) 

if no further constraints are imposed, and 

dM ~Af tot F 5 dFdA (1.6) 



Here and in the following we work in string units. 

5 After this paper was sent out, an interesting paper appeared |42| pointing out the possible existence 
of nonsupersymmetric AdS vacua at larger volumes, through a subtle balancing of nonperturbative and 
perturbative corrections in the potential |4ri| . 
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if one requires the cosmological constant A to be much smaller than F 2 . Scenarios 
in which supersymmetry breaking is driven by D-terms were also considered in 
|2U] , and it was pointed out that in the special case of supersymmetry breaking by 
an anti-D3 brane at the bottom of a conifold throat, low scales are more natural. 
Since we work in the large radius regime, there is no counterpart of this scenario in 
the M-theory compactifications we will study, so we will not get into details here. 
We should also point out that large classes of string compactifications have been 
proposed in j^S] where supersymmetry is broken at the KK scale. 

Next we summarise our results on the statistics of Freund- Rubin vacua. 

1.2 Freund-Rubin statistics 

Freund-Rubin vacua of M theory have geometry AdS 4 x X 7 , with X 7 a positively curved 
Einstein manifold, and can be understood as arising from the near-horizon geometry 
of A" coincident M5-branes, which become A" units of GVflux in the AdS-space. The 
compactification volume is fixed and depends on A" and the choice of X 7 . Typically, these 
geometries cannot really be considered as compactifications on X 7 in the usual sense, 
because the Kaluza-Klein scale tends to be of the same order as the AdS scale. 

Nevertheless we can study the distributions of AdS cosmological constants A and 
compactification volumes V. We will do this in section |21 for a model ensemble with 
X 7 = X/Zfe, where we vary k and N . This ensemble is extremely simple, and is therefore 
additionally useful as a simple toy model for testing counting methods. 

At fixed k ie for a fixed topology of the extra dimensions, we find the following distri- 
butions for A and V: 



This already shows a dramatic difference compared to IIB flux vacua. First, there are 
an infinite number of vacua, since A" is arbitrary. Secondly, the distribution of A is not 
uniform near zero, but diverges. Finally, large volumes are not suppressed — the larger 
N, the larger the volume becomes. 

Allowing both k and A" to vary ie by sampling the topology of the extra dimensions 
as well, these results significantly change. Now 



cW(A) ~ dA~ 2/ \ 



dM{V) ~ dV 6/7 . 



(1.7) 



dAf(A) ~ dA 



-2 



dJ\f(V) ~ dV 6 . 



(1.8) 
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The qualitative features are the same though: smaller cosmological constants and larger 
volumes are favored. 6 We also obtain the joint distribution for V and A: 



where 0(V,A) = 1 when V~ 3 < A < V~ 9/1 and zero otherwise (see also fig. HJ). One 
interesting feature that can be read off from this distribution is that at fixed A, V actually 
accumulates at smaller values, opposite to what we found for the unconstrained case. 
This is possible because the step function allows A to vary over a larger domain when 
V increases. This illustrates the importance of constraints for statements about which 
parameters are favored. Such issues become especially important if one wishes to interpret 
number distributions as probability distributions, since through such correlations, the 
dependence of these probabilities on one parameter may strongly influence the likelihood 
of values of the other. 

1.3 G 2 holonomy statistics 

The main part of this paper is devoted to the study of holonomy vacua. These have 

M = 1 supersymmetry and as usual come with many massless moduli. In 0j, a mechanism 

was proposed to stabilize all geometric moduli, basically by combining GVfluxes (which 

alone do not lead to any stable vacua 7 ) with a contribution to the potential induced by 

nonabelian degrees of freedom living on a codimension four singularity in the G<i manifold. 

In the low energy effective theory this has simply the effect of adding a complex constant 

c = c\+ic2 to the flux induced superpotential, given by a complex Chern-Simons invariant 

on the singular locus. The physics of this is presumably related to the Myers effect on 

D6-branes in type IIA, although this remains to be clarified. Anyway, regardless of its 

actual realization, this is also the simplest modification of the flux superpotential one can 

imagine and which leads to stable vacua. 

Unfortunately, no explicit constructions of G2 manifolds with the required properties 

are known (since G2 manifolds themselves are technically difficult to make), so we will 

just take these ensembles of flux superpotentials shifted with a constant as our starting 

6 Perhaps we should stress that by "favored" we do not mean "more probable". As emphasized earlier, 
we are computing number distributions at this level, not probability distributions. 

7 This is true in the classical supergravity approximation, in which we work throughout this paper. 

The results of ^] indicate that corrections to the Kahler potential may change this. 



(W(V, A) ~ 0(V,A) 



dV dA 
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point. Another problem is that not much is known about the metric on G2 moduli spaces, 
apart from the fact that they are derived from a Kahler potential K = —3 logV. To get 
around this problem, we follow two approaches. One, detailed in sectional is to keep the 
Kahler potential completley arbitrary and get general results. This will allow us to get 
quite far already. For some purposes however, more information is needed. In section HI 
we introduce an ensemble of model Kahler potentials which satisfy all known constraints 
for G2 moduli spaces. The models turn out to be exactly solvable, in the sense that all 
supersymmetric and nonupersymmetric vacua can be found explicitly, and we push the 
vacuum statistics analysis to the end. 

The results we find for general ensembles follows. 

1. The number of supersymmetric vacua in a region TZ of moduli space is, for large 02- 

Kusy ~4 3 Vol(ft). (1.10) 

Here 63 is the third Betti number of the G2 manifold, which is the number of moduli 
as well as the number of fluxes. Thus, vacua are distributed uniformly, at least in 
the large radius region of moduli space, where classical geometry and our analysis 
are valid. Note in particular that despite the absence of a tadpole constraint like in 
IIB, the number of vacua in the region of moduli space where our analysis is valid 
is finite. The result is quite similar to what was found for IIB flux vacua, with C2 
playing the role of the tadpole cutoff L. 

2. The distribution of volumes V (measured in lid Planck units) is given by 

dM = (kc2) b3 dV~ 3b3/7 . (1.11) 

Here k is a constant weakly decreasing with increasing 63. Thus, as in the IIB 
ensemble and its IIA mirror, large volumes are suppressed, and strongly so if 63 is 
large. An upper bound on V is given by 

V < (kc 2 ) 7/3 . (1.12) 

3. It can be shown that in a supersymmetric vacuum 




(1.13) 
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where A is expressed in 4d Planck units. The distribution of supersymmetric cos- 
mological constants is therefore entirely determined by the distribution of volumes: 



In particular for large 63, small cosmological constants are strongly suppressed, and 
bounded below by A > l/c^k 7 . This is radically different from the IIB case, where 
this distribution is uniform, and cosmological constants can be obtained as low as 
1/JVtot- The underlying reason for this difference is that in the IIB ensembles which 
have been studied, there are four times more fluxes than moduli, while in our G2 
ensemble there are only as many fluxes as moduli. This means there is much more 
discrete tunability in the IIB case, which makes it possible to tune the cosmological 
constant to a very small value. 

At this point the reader may wonder if this does not blatantly contradict what one 
would expect from string duality. This is not so, because our G 2 ensemble is not 
dual to the standard IIB ensemble with both NS and RR 3-form fluxes. On both 
sides, some of the flux degrees of freedom do not have a counterpart as flux degrees 
of freedom on the other side. In favorable circumstances these are dual instead 
to discrete geometric deformations away from special holonomy, although this has 
not been established in general [3B]- in any case, this shows that neither of the 
two ensembles is "complete". Conceivably, relaxing the G2 holonomy condition to 
something weaker could allow the cosmological constant in the M theory ensemble 
to be much more finely tuned, and one could imagine that one would then get a 
uniform distribution again, although this is by no means certain. It would be very 
interesting to verify this, but at present a detailed description of such would be 
vacua is unknown. 

4. The distribution of supersymmetry breaking scales is similarly quite different from 
the IIB case, again essentially because of the lack of discrete tunability. Because 
there are as many fluxes as moduli, one can express all flux quanta as a function 
of the susy breaking parameters at a given point in moduli space. The equation 
V — at that point thus becomes a complicated quadratic equation in the F if 
whose solutions, apart from the supersymmetric F = 0, are not tunably small. This 
means the analysis of j20] does not apply, and general, complete computations of 
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distributions become hard. 

A number of simple observation can be made though. Since we have a system of 63 
quadratic equations in 63 variables F iy the number of nonsupersymmetric branches, 
for a fixed choice of flux can be up to 2 b:i (minus the supersymmetric vacuum). 
In the model ensembles of section 01 we show that this number can actually be 
obtained. Since essentially all of these break supersymmetry at a high scale, this 
is further support for the idea that string theory has many more F-breaking flux 
vacua with high scale breaking than with low scale. 

More precisely, we show that in these ensembles, any vacuum with cosmological 
constant A ~ (assuming it exists) has supersymmetry breaking scale 

M L y = ^rn 2 p , (1.15) 

where r « 0.1. Thus, as for the cosmological constant, the scale of susy breaking 
is set by the volume V, and similarly its distribution will largely be determined by 
the distribution of V. Using the volume bound ()1.12j) . 8 this gives a lower bound on 
the supersymmetry breaking scale: 

M s 2 usy > Ar^H (1-16) 

For moderate values of c 2 (recalling that k is decreasing with increasing 63), the 
supersymmetry breaking scale will therefore always be high in these ensembles. 

This behavior is quite different from Type IIB and the generic ensembles studied 
in [20], which again can be traced back to the fact that there is limited discrete 
tunability because there are only as many fluxes as moduli. It is conceivable that 
adding discrete degrees of freedom, by allowing deformations away from special 
holonomy, would bring the distribution of supersymmetry breaking scales closer to 
the results of [20] • in any case this would not change the qualitative conclusion that 
higher superymmetry breaking scales are favored. 

Finally, in section|3]we present and study a class of models defined by Kahler potentials 

which give a direct product metric on moduli space. Though very rich, the semi-classical, 

8 which was strictly speaking derived for supersymmetric vacua, but the distributions of supersym- 
metric and nonsupersymmetric vacua over moduli space are expected to be very similar, as we confirm 
in the model ensembles of section 0] 
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supergravity vacuum structure of these models can be solved exactly, allowing us to 
explicitly verify our more general results. Vacua are labeled by (A^,o"j), i = 1, ...,63, 
where the N{ are the fluxes and o~j = ±1. Putting all cr.j = +1 gives a supersymmetric AdS 
vacuum, the 2 b:i — 1 other choices correspond to nonsuperymmetric vacua. As the fluxes are 
varied, we find that these are uniformly distributed over the moduli space and that there 
are roughly an equal number of de Sitter vs. anti de Sitter vacua. Not all of these vacua 
exist within the supergravity approximation and we analyse the conditions under which 
they do, finding that an exponentially large number survive. Finally, after analysing the 
stability of these vacua we found that all de Sitter vacua are classically unstable whilst an 
exponentially large number of non-supersymmetric anti de Sitter vacua are metastable. 
As expected on general grounds, the supersymmetry breaking scale is typically high. 

These model ensembles are reminiscent of the effective field theory landscapes recently 
considered in |2E|. In particular, at large 63, the distributions we find are sharply peaked, 
see e.g. fig. |H]and fig. for examples of distributions of cosmological constants and super- 
symmetry breaking scales. The cosmological constant in these ensembles does not scan 
near zero; there is a cutoff at A ~ — c^ 5 . It is conceivable though that this would change 
for more complicated Kahler potentials which do not lead to direct product metrics. 



In this section we will study the statistics of Freund-Rubin vacua. We will study a model 
ensemble for which the basic distributions of cosmological constants and volumes are 
extremely simple compared to say IIB flux vacua or GVholonomy vacua. After a brief 
review of Freund-Rubin vacua, included for completeness, we describe the distributions. 
Note that, since 4d Freund-Rubin vacua are dual to three dimensional CFT's, we are 
equivalently studing the statistics of such theories. 

2.1 Freund-Rubin basics 

Freund-Rubin vacua [33] are near horizon geometries of brane metrics in which there is 
an AdS factor. We will focus on the M2-brane case for which the Freund-Rubin metric 
takes the form 



2 Freund-Rubin Statistics 




(2.1) 
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Here AdS a is four dimensional anti de Sitter spacetime with radius a and X is a compact 
7- manifold, and go(X) is an Einstein metric on X with positive scalar curvature.: 



Since the actual metric on the seven compact dimensions is rescaled by a 2 , the scalar 
curvature of X and the four dimensional AdS is of order a -2 . 

This means that if the volume of X as measured by go(X) is order one compared 
to that of a round unit 7-sphere, then the masses of Kaluza-Klein modes are of order 
the gravitational mass in AdS (i.e. the inverse AdS radius). In this case, there is no 
meaningful low energy limit in which one can ignore the dynamics of the Kaluza-Klein 
particles. 

However, if the radii of X in the go{X) metric are smaller than one [SU], there will 
be a mass gap between the gravitational fluctuations and Kaluza-Klein modes. In this 
situation there is a meaningful low energy limit. 

Freund-Rubin vacua can be both supersymmetric or non-supersymmetric. In fact, 
for every M = 1 Freund-Rubin vacuum there exists a semi-classically stable, non-super- 
symmetric Freund-Rubin vacuum. This is obtained by simply changing the M2-branes 
to anti M2-branes, or equivalently by changing the sign of G. This was known in the 
Kaluza-Klein supergravity literature as "skew- whiffing" In fact, there are presumably 
many more non-supersymmetric Freund-Rubin vacua, since the Einstein equations with 
positive curvature are a well posed problem, and this is irrespective of super symmetry. 
Unfortunately, for a general, non-supersymmetric Freund-Rubin solution, semi-classical 
stability of the vacuum can only be determined on a case by case basis. The "skew- 
whiffed" solutions are always stable. 

Finally, it was pointed out only recently that Freund-Rubin vacua in four dimensions 
can have chiral fermions |5U]. The mechanism for this — which evades Witten's no go 
theorem [51] is that the chiral fermions are wrapped membranes localised at singularities, 
precisely as in the G^-holonomy case studied in EH]. The possibilities for realistic 
particle physics phenomena in this context are little explored at the present time and 
certainly deserve further investigation. Some suggestions for the construction of more 
realistic models were given in [50J. 




(2.2) 
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2.1.1 Quantization of a 

The Freund-Rubin metric has one parameter a. On the other hand, brane metrics also 
have one parameter, N the number of branes. Therefore, the two parameters are related. 
The relation is 

a" = £ (2.3) 

where Vq is the number of Planck volumes of X as measured by go(X). A Planck volume 
is L 7 p where L p is the eleven dimensional Planck length. 

To show that (|2.3|) is correct we need the formula for the G-flux in the vacuum. This 

is 

G~ -rfVol(AdS a ). (2.4) 
a 

The charge or number of branes is measured by integrating the flux *G around the brane, 
i.e. over X: 

N = [ *G = a 6 V . (2.5) 
Jx 

This reproduces ()2.3|) . We will write physical quantities in terms of N rather than a. 
The four dimensional Planck scale m p is related to the eleven dimensional M p via 

2 7.T 7 

-E- „ aV ~ (2.6) 

M p V* 

This is the volume of X in lid Planck units. The vacuum energy is 



M 1 V 



A £mj ~ -^ml (2.7) 



2.2 Statistics 



Obviously at the classical level is an arbitrary integer and there are hence an infinite 
number of vacua. Clearly however as A^ tends to infinity the volume of X and the radius 
of AdS also tend to infinity, and the cosmological constant goes to zero. A more sensible 
question is therefore how many vacua are there with compactification volume Vx < V^? 
Since Vx ~ ^t/^, V* is ^4jf- So this question has the simple answer: 

KJYx < V.) = iV max ~ V l/7 V^ 7 . (2.8) 
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k 



V 



Figure 1: Left: Lattice of vacua in (k, iV)-space. The green solid lines have constant V and the red 
dashed lines constant A -1 . Both are increasing with N and k. Right: Vacua mapped to (V, A _1 )-space 
with V < 200, A^ 1 < 2000. The lower and upper boundaries correspond to k = 1 resp. N = 1. 

Similarly the number of vacua with |A| < A* is infinite. However, since |A| ~ Vq /2 N- 3 / 2 
the number of vacua with |A| > A* in 4d Planck units is 

Kc(\M > A*) ~ V 1/3 Af/\ (2.9) 

These results are valid for a fixed Einstein manifold with normalized volume Vq. We 
can also ask: how do these numbers change if we vary the Einstein manifold XI We 
can construct simple examples of ensembles of vacua in which X itself is also varying as 
follows. 

If X has a U(l) symmetry then we can quotient by a discrete subgroup C U(l). The 
quotient X/Z k has volume ^o(^) = , but is still Einstein with the same cosmological 
constant. There are many examples of such families. A very simple one has X = S the 
round 7-sphere, regarded as the sphere in C 4 . This solves the equations of motion. If 
the coordinates of C 4 are denoted by Zj, then : Zi — > e^Zi acts freely on S 7 . The 
quotients |- provide a family of topologically distinct Freund-Rubin solutions labeled by 
k. 
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For such a set of vacua labeled by (N, k) we have (from the formulas in section E.l.lj) 9 

kN 3 = A" 2 (2.10) 
kN 7 = V 6 . (2.11) 

We may now ask what the number of vacua is in a given region of (V, A)-space. This 
amounts to counting lattice points in the corresponding region in (k, A r )-space, as illus- 
trated in fig. H For sufficiently large regions, an estimate of this number is given by the 
area in (k, AQ-space. Such an estimate is in general better for higher dimensional lattices 
and convex regions, but gives a reasonable approximation in this two-dimensional case as 
well, at least for sufficiently large regions, as we will demonstrate below. 

Instead of computing the area in (k, N) space, one can integrate a suitable vacuum 
number density over the region in (V, A)-space. This density is given by 

dU = dkdNQ(k — 1) Q(N — 1) 

= d(v^ 2 A- 7 / 2 ) d(y 3/2 A 1/2 ) e(^ 9 / 2 A- 7 / 2 - 1) e(y 3 / 2 A 1 / 2 - 1) 

° dVdAQ(l-V 9 A 7 )Q(V 3 A-l) (2.12) 



V A A 4 
3(A' X ) 2 



dVdA- 1 6(A- 1 - V 9/7 ) 0(V 3 - A' 1 ) (2.13) 



V 4 

Here @(x) denotes the step function, equal to 1 if x > and otherwise. The last line 
is included for comparison with the discrete density in the (V, A -1 )-plane shown in the 
figure. 

We see for example that at fixed A, the density of vacua dM ~ dV/V 4 is higher 
towards smaller volumes. However, this does not mean that in total there are more vacua 
at smaller volumes, as the step functions cut out an allowed region for A that grows when 
V increases. Going back to the original parametrization in terms of k and A^ (which is 
more useful to compare to the exact discrete distribution), we get that in the continuum 
approximation, the total number of vacua with bounded volume is given by 

/oo i>oo 
dN J dkQ(j^-k) (2.14) 

where the step function enforces the condition that we count A r 's and k's such that the 

volume is bounded by V* (i.e. the area under a solid green line in fig. Q). This integral 

9 We are ignoring numerical coefficients which are irrelevant for our applications, and will abuse nota- 
tion and write A instead of |A| in what follows. 
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gives the answer 

V 6 7V 6/7 

K ac (v x < K) = -t + 1 - —r- ~ K 6 - (2.15) 

o o 

The subleading terms can be dropped since the continuum approximation requires V* 3> 1. 
The leading term here arises from the small N region so most of the vacua are here, as 
can be seen in the figure as well. So we find that the total number of vacua actually grows 
with the volume, in contrast to the decreasing dV/V A density at fixed A. Note that this 
V® growth is also quite different from the number of vacua with Vx < V* at fixed k, given 
by (|2.8p . which scales as V^ 7 . Therefore, sampling topologies of the extra dimensions 
significantly changes the distribution of vacua. 

These observations illustrate the obvious but important fact that the distribution of 
a quantity can depend strongly on the constraints and ensembles considered. 

Another important issue is the validity of the continuum approximation. Thanks to 
the simplicity of this ensemble, we can compare our result (|2.15|) with the exact answer: 

NtJYx < v.) = - k ) « E w ( 2 - 16 ) 

fc,iV>l 1<N<V? /7 

The approximation we made here (neglecting integrality of k max ) becomes asymptotically 
exact in the large V* limit. In the same limit, the sum over N gives a factor C(7), so 

N vac (v x < v;)^c(7)K 6 - (2-17) 

Since £(7) ~ 1.008, this is about a factor 6 larger than the continuum result (J2.15)) . The 
reason for this discrepancy is the fact that the region under consideration in the (k, N)- 
plane is very elongated along the fc-direction, and most contributing lattice points are at 
the boundary N = 1 (see fig. HJ, so the integral is not a very accurate approximation of 
the sum. The continuum approximation will be better for regions in (V, A)-space that 
stay away from the boundary. Even in the case at hand though, the power of V* itself is 
correct, and since the coefficient is irrelevant at the level of accuracy we are working at 
anyway, the continuum approximation is satisfactory for our purposes. 

Similarly we could ask how many vacua are there with A > A* (the other way around 
obviously gives a divergent answer). This is 

/oo poo 1 

dNj dke(j^-k)^A^ (2.18) 

Again this is a stronger growth than the fixed k distribution (|2.9j) . 
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^From the joint distribution (j2.12|) . one gets the distribution of one variable given the 
other one. As noted before, at fixed A, vacua accumulate towards the lower bound on 
the volume V, which is opposite to what one has for the distribution without constraint 
on A. For instance, when A~ 2 < V® and A* 2 ^ 3 > V^ 1 , one finds that Af vac (V > V*; A > 
A*) ~ A^T 3 ^ 3 which favors both small volume and cosmological constant. 

In summary: we have shown that sampling the topology of the extra dimensions dras- 
tically changes the distribution of cosmological constants and volumes (i.e. eleven di- 
mensional Planck scales). A priori, Freund- Rubin vacua statistically congregate in large 
volume, small cosmological constant regions. However, we saw that strong correlations 
between cosmological constant and volume exist, causing the effective distribution of one 
quantity to depend significantly on the constraints imposed on the other. For example 
at fixed cosmological constant, vacua actually prefer the lower volumes within the al- 
lowed range. Finally, we stressed that there exist at least as many non-supersymmetric, 
meta-stable Freund-Rubin vacua than supersymmetric. 

One important point to keep in mind is that we have sampled a very simple set of 
Einstein manifolds. One would certainly like to know how representative the distribution 
of volumes and vacuum energies for the (N, k) ensembles we studied is in the ensemble of 
all Einstein manifolds with positive scalar curvature. Some more complicated ensembles 
of Einstein 7-manifolds are described in loTi] and jSS] gives a formula for the volumes 
which could be used to study the distributions of vacua as we have done here. 

3 G 2 Holonomy Statistics I: general results 

G2-holonomy vacua are compactifications of M theory to four dimensions which, in the 
absence of flux classically give 4d M = 1 vacua with zero cosmological constant. These 
classical vacua have 63(A) complex moduli, of which the real parts are axions t« and the 
other half Sj are the massless fluctuations of the metric on X. 

The addition of fluxes when X is smooth does not stabilise these moduli, as the induced 
potential is positive definite and runs down to zero at infinite volume. However, if X has 
an orbifold singularity along a 3-manifold Q, additional non-Abelian degrees of freedom 
arise from massless membranes j^Zj. Nonabelian flux for these degrees of freedom then 
gives an additional contribution to the potential which can stabilise all the moduli if Q 
admits a complex, non-real Chern-Simons invariant This is the case if, for instance, 
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Q is a hyperbolic manifold. 

The vacua studied in [3] were supersymmetric with negative cosmo logical constant. 
In fact, it was shown that in the large radius approximation, for a given flux within a 
certain range, there is a single supersymmetric vacuum (in addition to an unstable de 
Sitter vacuum). In principle however there could be other, non-supersymmetric, vacua 
and one of our aims here is to study this possibility. One might wonder if any of these 
vacua could be metastable de Sitter. We will answer this question to a certain extent. 

One of the difficulties in studying GVholonomy compactifications is that G^-holonomy 
manifolds are technically very difficult to produce. For instance, we still do not know 
whether or not there exists a GVholonomy manifold with a non-real Chern-Simons in- 
variant. So how can we hope to study the statistics of such vacua? As we will see below 
the superpotential of these GVcompactifications with flux is very simple and does not 
contain much information about X. Instead, this information comes through the Kahler 
potential on moduli space, which could be a quite complicated function in general, of 
which very little is concretely known. One approach then is to try to obtain general re- 
sults for an arbitrary Kahler potential. We will do this in the next section by extending 
some of the general techniques which were developed in the context mainly of IIB flux 
vacua in EJ 120] • Secondly we could study particular ensembles of model Kahler poten- 
tials and hope that the results are representative in general. We will follow this approach 
in section where we study a particular class of model Kahler potentials which allow 
explicit construction of all supersymmetric and nonsupersymmetric vacua. 

3.1 G 2 compactifications with fluxes and Chern-Simons invariants 

Let X be the G2 holonomy compactification manifold. The complexified moduli space M. 
of X has dimension n = b 3 (X) and has holomorphic coordinates z l , defined by 




(3.1) 



20 



where (p is the (^-invariant 3-form on X and £m is given by l/2K,f 1 = 2tc/£ 9 m . 10 The 
metric on M. is Kahler, derived from the Kahler potential 11 [51] 

K{z,z) = -3\n(^V x {s)), (3.2) 

where Vx = Vol(X)/£ 7 M is a homogeneous function of the s l of degree 7/3. This classical 
metric will receive quantum corrections, but at large enough volumes such corrections can 
be argued to be small. 
We turn on 4-form flux 

G/t M = N lP \ (3.3) 

where JVj £ Z and p % is a basis of if 4 (X, Z) , and also assume the presence of a complex 
Chern-Simons contribution as described above. This induces a superpotential (551 1551 I5U 

El 

W(z) = ^(N i z i + ci + ic 2 ) (3.4) 

on Ai, where C\ and C2 are the real and imaginary parts of the Chern-Simons invariant. 

The corresponding potential is obtained from the standard four dimensional super- 
gravity expression (for dimensionless scalars): 

V = 4e K [fFiFj - 3\W\ 2 ^ , (3.5) 

where 

Fi = D z iW = (d z i + d z ,K)W = Ni + ^d s ,KW. (3.6) 

We will put k 4 = 1 in what follows. 

Since we are working in the large radius regime, the axions t l essentially decouple 
from the moduli s l , and all nontrivial structure resides in the latter sector. This is seen 
as follows. Writing 

W = W 1 (t)+i W 2 (s), Ki = d si K(s), Kij = d si d si K = Ag fj (3.7) 
and so on, the potential (|3.5|) becomes 

V = e K (4K ij (Re Fi)(Re Fj) + K^KiKjW? -3Wi) 

= e K (4K ij (Re F;)(Re F 5 ) - 3W 2 (s) 2 + AW^t) 2 ) . (3.8) 

10 In terms of £m, we have T2 = 2ir/£\j, T5 = 2ir/£^, and / G/£\j € Z. Thus, the instanton action is 

e 2 ™ 2 and t has periodicity 1. 

11 Our normalization conventions for z, W and K are slightly different from |54l . The value of the 

normalization coefficient appearing in (|3.2(l in front of Vx is verified in Appendix A. 
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In the last line we used the fact that the volume Vx is homogeneous of degree 7/3, which 
implies the following identities: 

KiJ = -7, K iiS ? = -Ki. (3.9) 

The second is obtained from the first by differentiation. 

Since W± = Nit 1 + c\ and everything else in (|3.8|) depends only on s, it is clear that 
any critical point of V will fix 

N% + c 1 = Q (3.10) 

and therefore W\ = Im F — 0. Apart from this, the f are left undetermined, and they 
decouple from the s\ From now on we will work on this slice of moduli space, so we can 
write 

V = e K (4K ij FiFj - 3W 2 2 ) (3.11) 

with 

Fi = D si W 2 = (d s i + X -K t )W 2 = N, + ^K t (N jS > + c 2 ). (3.12) 

The geometry of the real moduli space Ai parametrized by the s l is the real analog of 
Kahler, often called Hessian, with metric Kij. 

3.2 Distribution of supersymmetric vacua over moduli space 

The aim of this section is to find the distribution of vacua over Ai, along the lines of 
[HI El 120] • The condition for a supersymmetric vacuum in the above notations is 

D s iW 2 (s) = (3.13) 

In what follows we will drop the index 2 to avoid cluttering. The number of solutions in 
a region 1Z of Ai , for all possible fluxes G but at fixed c 2 , is given by 

N snsj (c 2 ,TZ)= [ drs^DiWydetdiDjWl (3.14) 

The determinant factor ensures every zero of the delta-function argument is counted 
with weight 1. An approximate expression A4usy f° r the exact number of vacua iV susy is 
obtained by replacing the discrete sum over iV by a continuous integral, which is a good 
approximation if the number of contributing lattice points is large (which will be the case 
for sufficiently large c 2 ). Thus 

A4us y = J d n s J d n N6 n {D i W)\detd i D j W\. (3.15) 

22 



Differentiating (j3.12|) and using (j3.13|) and ()3.9|) . one gets 

diDjW = (^Kij - \k { K^ W = l -K ik [b] + \s k K^j W (3.16) 

One eigenvector of the matrix in brackets is s\ with eigenvalue —5/2 (this follows again 
from (J3.9|0 . On the orthogonal complement {v\Kjv^ = 0}, the matrix is just the identity, 
so all other eigenvalues are 1. Thus, 



det d t D,W\ = det K l3 (3.17) 



Furthermore, by contracting ()3.12|) with s\ we get at Fi = 0: 

W = -\c 2 . (3.18) 

To compute 1)3.15)1 . we change variables from iVj to F{. The Jacobian is 

Jac = | detd^FA- 1 = | det(« + ^K.s 1 )^ 1 = 2/5. (3.19) 

J 2 

Putting everything together, we get 

Kusy = J d"s J d n F5 n (F)(^j detKtj (3.20) 

d n s det K (3.21) 



C 2 



5 



/i \ n 



'/v 



Vol(ft) (3.22) 

where 7?. is the part of the complexified moduli space projecting onto 71 (i.e. the direct 
product of 1Z with the n-torus [0,1]™ swept out by the axions f), and the volume is 
measured using the Kahler metric gq = Kij/A on this space. 

Thus, supersymmetric vacua are distributed uniformly over moduli space. This result 
is similar to the Type IIB orientifold case studied in pj, but simpler: the Type IIB 
vacuum number density involves curvature terms as well, and in order to get a closed 
form expression for n > 1, it was necessary there to count vacua with signs rather than 
their absolute number. 

Note that in any finite region of moduli space, the number of vacua is finite, because 
C2 is finite despite the absence of a tadpole cutoff on the fluxes. In particular, the total 
number of vacua in the large radius region of moduli space (where our computation is 
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valid) is finite. For IIB vacua on the other hand, finiteness is only obtained after imposing 
the tadpole cutoff J F A H < L* = x/24 on the fluxes. In a way, the Chern-Simons 
invariant c 2 plays the role of here. 

Example 

The simplest example is the case n = 1. Then homogeneity fixes Vx ~ s 7//3 , so K ss = 
and 

A/" susy (c 2 , s > s*) = ^ / ds — = -—. (3.23) 

Thus, vacua become denser towards smaller volume of X, and from this equation one 
would estimate there are no vacua with s > 7c 2 /5, which is where Af susy drops below 1. 
Indeed, it is easily verified that the exact critical point solution is s = —7c 2 /5N, so the 
largest possible value of s, obtained at iV = —1, is precisely 7c2/5. Better even, from the 
explicit solution it easily follows that the approximate distribution (J3.23)) becomes in fact 
exact by rounding off the right hand side to the nearest smaller integer. 



3.3 Large volume suppression 

Although the precise form of the metric for n > 1 is unfortunately not known, one general 
feature is easy to deduce: flux vacua with large compactification volume are suppressed, 
and strongly so if the number of moduli is large. 12 This follows from simple scaling. 
If s — > As, K shifts with a constant, so Kij — > \~ 2 Kij and the measure d n sdetK — > 
A~ n d n s det K. Thus we have for example 

Mu Sy (c 2 , j > S*) = (C 2 /S*)" A4usy(l, > 1), (3.24) 

and 

A4us y (c 2 , V x > K) = K (c 2 /K 3/ T , (3.25) 

12 This observation was originally made in collaboration with M. Douglas in the context of (mirror) IIB 
flux vacua |44| . 
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and 



where k n is independent of c 2 and K- 13 Hence large compactification volume is strongly 
suppressed when n is large. 

To get an estimate for the absolute number of vacua with Vx > K, one would need to 
estimate k n in 1)3.25)1 . This would be quite hard in general even if the metric was explicitly 
known, so to get an idea let us do this for a simple toy model, taking the volume to be 
the following homogeneous function of degree 7/3: 

M«) =(£«?) • (3-26) 

Then K = -|log(s 2 ), 

^ = ^(2^-^), (3.27) 

\detK\ = (J^J . (3.28) 

This metric is actually not positive definite for n > 1 and therefore unphysical, but let us 
proceed anyway. A more sensible model will be presented in the next section. We have 
Vx > V* iff \s\ > s*, with = V^ 7 . Thus 

Ku Sy (c 2 ,V x >V*) = f d n s^ (3.29) 

(3.30) 

where Q n -i = r(n/2) * s ^ ne area °^ ^ ne ( n ~~ l)-dimensional sphere of radius 1. So 

M m ,M,V X >V.) = (fj w - (3.31) 



7\/2ire C2 V / 6 c 2 



(3.32) 



m 

j 



5 s*y/n ) \\fns 

In going from the first to the second line we used Stirling's approximation m! ~ (m/e) 
valid for large m. u 

13 Even though the obvious metric divergence at Vx — is avoided by bounding Vx from below, it 
might still be possible that this bound alone does not determine a finite volume region in s-space. Then 
the left hand side will be infinite and the scaling becomes meaningless. In this case, additional cutoffs 
should be imposed, which will complicate the dependence on V*, but large volume suppression is still to 
be expected. 

14 In fact ml > (m/e) m , so the above approximation for Sf BWy gives an upper bound. 
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This confirms the general result based on scaling, but now we also have information 
about the absolute numbers. In particular, this formula suggests all vacua for this toy 
model must have compactification volume approximately bounded by 



Similar bounds can be expected for other models, as we confirm in section 0] In general, 
large extra dimension scenarios with, say, micrometer scale compactification radii are 
therefore excluded in these ensembles unless C2 is exceedingly large. 

3.4 Distribution of supersymmetric cosmological constants 

The vacuum energy in a supersymmetric vacuum is, in 4d Planck units, using ([3. 18)) : 



where a = 3/400 tc « 0.002. Therefore, the only way to get a small cosmological constant 
is to have large compactification volume. This contrasts with the ensemble of type IIB 
flux vacua, where very small cosmological constants can be obtained at arbitrary complex 
structure (i.e. arbitrary mirror IIA volume). The underlying reason is the fact that in IIB, 
there are four times as many fluxes as equations D$V = 0, so at a given point, there is 
still a whole space of (real) fluxes solving the equations. This freedom can be used to tune 
the cosmological constant to a small value. Here on the other hand, there is only one flux 
per equation, so at a given s, the fluxes are completely fixed, and no freedom remains to 
tune the cosmological constant. This would likely change however if more discrete data 
were turned on, such as the M theory duals of IIA RR 2-form flux or IIB NS 3-form flux. 
Unfortunately these are difficult to describe systematically in M theory in a way suitable 
for statistical analysis. 

Let us compute the distribution of cosmological constants more precisely. Equation 
(I3.34J) implies that this follows directly from the distribution of volumes. Using (|3.25j) . 
we get 




(3.33) 




x 



(3.34) 




(3.35) 



The corresponding distribution density is therefore 



dM/d\ ~ A (n - 7)/7 . 



(3.36) 
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In particular, for n < 7, the distribution diverges at A = 0, while for n > 7, the density 
goes to zero. In the toy model, from ()3.33|) . we get an expected bound: 



IAI > W - (3.37) 



For large n, this is much larger than the naive l/J\f susy which was found to be a good 
estimate in the Type IIB case, where the cosmological constants of supersymmetric vacua 
are always distributed uniformly near zero 

Because Gi manifolds with many moduli are much more numerous than those with 
only a few, we can thus conclude that small cosmological constants are (without further 
constraints) strongly suppressed in the ensemble of all supersymmetric G2 flux vacua. 

3.5 Nonsupersymmetric vacua 

A vacuum satisfies V = 0, and metastability requires V" > 0. Thus, the number of all 
metastable vacua in a region 1Z is given by 



N vac = Y) / d n s 8{V) I det V"\ e(V"). (3.38) 



In principle one could again approximate the sum over iV by an integral and try to solve 
the integral by changing to appropriate variables, as we did for the supersymmetric case 
in section El and as was done in [20J for supersymmetry breaking scales well below the 
fundamental scale. In practice, we will encounter some difficulties doing this for G2 vacua. 
By differentiating equation (|3.11j) . one gets 

diV = e K (8(D i D j W)D j W - 6 WD t W) (3.39) 

where Di denotes the Hessian Levi-Civita plus Kahler covariant derivative. The matrix 
Mij = DiDjW is related to the fermionic mass matrix, and with this notation the critical 
point condition becomes 

Mij F j = ^WFi. (3.40) 

In j2ni 5 a similar equation was interpreted at a given point in moduli space as a linear 
eigenvalue equation for the supersymmetry breaking parameters F, assuming the matrix 
M (denoted Z there) to be independent of F. This is indeed the case for the Type IIB 
flux ensemble, but it is not true for the G2 ensemble. The reason is again the fact that 
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there are only as many fluxes as moduli here, so a complete set of variables parametrizing 
the fluxes Ni is already given by the Fi (affmely related to the Ni as expressed in (|3.12Jl ). 
Therefore, all other quantities such as W and and so on must be determined by the 
Fi. Indeed, a short computation gives: 

W = ~l(s l F t + c 2 ) (3.41) 
5 

Mij = ^((Kij — ^KiKj)W + KiFj + KjFi — K^F k ) (3.42) 

This means that at a given point in moduli space, ()3.40j) is actually a complicated system 
of quadratic equations in F, with in general only one obvious solution, the supersymmetric 
one, F — 0. In total one can expect up to 2 n solutions for F. The nonsupersymmetric 
solutions will generically be of order C2. 15 In particular it is not possible to tune fluxes 
to make F parametrically small, so it is not possible to use the analysis of [20] here, and 
there is no obvious perturbation scheme to compute ()3.38|) . 

Example 

As a simple example we take the case n — 1. As noted in section l3~2*l homogeneity 
then forces K = —71ns, hence M ss = (7c 2 — 5Fs)/2s 2 , and ()3.4(Jj) (considered as an 
equation for F at fixed s) has solutions F = and F = ^j 2 -- Neglecting the metastability 
condition, the continuum approximated number density of nonsupersymmetric vacua is 
given by 

f 14 C9 14 Co 

Nonsusy = / dN 5(F - —±) \d s (F - —^)\ds 

= f 2 -dFS(F- S ^&s S 
J 5 K s ' 2s 2 

= ^ds. (3.43) 

s 2 

We changed variables from iV to F in the integral using the Jacobian (J3.19)) . This result 
should be compared to the supersymmetric density implied by ()3.23j) : 

dA/nonsusy = 5 dAf snsy . (3.44) 

This does not mean that for a given flux, there are on average five nonsupersymmetric 
critical points for every supersymmetric one. In fact, as we will see below, and as was 
15 The physical supersymmetry breaking scale has an additional factor e K / 2 . 
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already pointed out in for n — 1, there is exactly one nonsupersymmetric critical point 
for every supersymmetric one: the supersymmetric minimum is separated from s = oo by 
a barrier, whose maximum is the (de Sitter, unstable) nonupersymmetric critical point. 
Equation ()3.44j) only expresses that in a given region of moduli space (at large s % ), there 
are on average 5 times more nonsupersymmetric vacua then supersymmetric ones. This is 
simply because in the region of moduli space under consideration, the nonsupersymmetric 
critical points are located at five times the value of s of the supersymmetric critical points. 



3.6 Supersymmetry breaking scales 

For more moduli, things become more complicated. A few useful general observations 
can be made though. If we require V = (which remains a good approximation for what 
follows as long as \V\ <C e K c\) we can make a fairly strong statement about the value 
of the supersymmetry breaking scale. Contracting (|3.40|) with s % and using (|3.41j) - ([3.42|) 
and (|3.9j) . we get 

25 F 2 - 3(s • F) 2 — 8 c 2 (s • F) = 0. (3.45) 

Together with V ~ 4F 2 — 3W 2 = 0, this gives a system of two equations in two variables, 
s ■ F and F 2 , with solutions: 

F 2 = M s-F= 3 f. (3.46) 
The physical supersymmetry breaking scale for such vacua (assuming they exist) is 



9 9 / — ~ — ~. a/3 c 2 m 2 Mlc 2 

ML y - m^AeKKvF^ = — — | ~ -f-, (3.47) 



where m p = 1/k 4 and M p are the 4- and 11-dimensional Planck scales. 

Let us plug in some numbers to get an idea of the implications. If we identify M v 
with the unification scale M un if ~ 10 16 Gev, this means that M susy ~ ^/cq. 10 14 ' 5 GeV, and 
the gravitino mass M g ~ M 2 usy /m p ~ c 2 10 10 GeV. Since c 2 is at least of order 1, this 
estimate implies (under the given assumptions) that in this ensemble supersymmetry is 
always broken at a scale much higher than what would be required to get the electroweak 
scale M cw ~ 100 GeV without fine tuning. 

In fact a slight extension of this calculation shows that even if one allows the addition 
of an arbitrary constant to the potential to reach V = (e.g. to model D-terms or 
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contributions from loop corrections), the supersymmetry breaking scale in this ensemble 
is still bounded from below by the scale C2M 3 /m p . 

One could of course question the identification M p ~ M uni f. Lowering the lid Planck 
scale down to M p ~ 10 13 Gev for example (which requires Vx ~ 10 12 ) gives a supersym- 
metry breaking scale M susy ~ y/c^lO 10 Gev and gravitino mass m g ~ C2 10 Gev, which for 
C2 not too large would give low energy supersymmetry. 

Note however that in analogy with the supersymmetric case, and based on general 
considerations, we expect suppression of vacua at large volume and therefore also sup- 
pression of low energy supersymmetry breaking scales in this ensemble. This is further 
confirmed by the exactly solvable models we will present in section 0] In particular we 
get a lower bound on the supersymmetry breaking scale from the upper bound on the 
volume Vx, which in general we expect to be of the form (|3.33j) . i.e. Vx < (c2/r n ) 7 ^ 3 , with 
r n weakly growing with n. Using the relation (J3.47|) between M susy and Vx, this implies 
the following lower bound for the supersymmetry breaking scale 

7/4 

M snsy /m p > (3.48) 

C 2 

Getting M susy below 10 12 Gev in the case of many moduli would thus require C2 to be at 
least of order 10 6 . 

Using the volume distribution (J3.25|) . we furthermore get an estimate for the distribu- 
tion of supersymmetry breaking scales (in 4d Planck units): 

cW„y ~ h c 2 n/7 d(M*:J y 7 )/m^ 7 , (3.49) 

where k n is a constant independent of C2 and M*. 

Here we have not yet taken into account the tuning required to get a tiny cosmological 
constant: |A| ~ |4F 2 — 3iy 2 |/V|- < A*. For a given volume Vx (or equivalently a given 
supersymmetry breaking scale), this requires tuning the fluxes such that |4F 2 — 3W 2 \ < 
A*^, which can be expected to at least add another suppression factor A^V^ ~ A*/M s 4 usy . 
The suppression may in fact be stronger, if the distribution of cosmological constants is 
not uniform but more like a Gaussian sharply peaked away from A = 0. Such distributions 
are quite plausible in these ensembles, as will be illustrated by the model ensemble we 
will study in section 0] Another potentially important factor which we are neglecting in 
this analysis is the metastability constraint (this was found in [20] to add another factor 
M s 4 to the distribution in the ensembles studied there). 



30 



Finally, when one also takes into account the observed value of the electroweak scale 
M ew , there is an additional expected tuning factor presumably of order Mg W mp/M s 4 usy (in 
the region of parameter space where this is less than 1) ^TJEHEl- Putting everything 
together, this gives (for > M s 2 usy > M ew m p ): 

dN~ k n KMlcf /7 \dMS 7 -^\/m^ 7 - 2 (3.50) 

So we see that for n < 14, the Higgs mass and cosmological constant tunings tilt the 
balance to lower scales, while for n > 14 higher scales are favored, and strongly so if n is 
large. For n < 14 we should keep in mind however that there is an absolute lower bound 
on the supersymmetry breaking scale, given by (|3.48|) . which will further be increased 
by the additional tunings of cosmological constant and Higgs mass. We should also not 
forget that this is only a naive analysis; in principle a full computation of the measure 
should be done along the lines of (201 , but as we discussed, this does not seem possible in 
the present context, because of the absence of a small parameter. 

Nevertheless, the above consideration indicate clearly that low energy supersymmetry 
is typically disfavored in G2 flux ensembles, and even excluded if C2 is less than 10 6 . 

4 G 2 statistics II: model Kahler potentials and exact so- 
lutions 

In the previous section we obtained a number of general results about distributions of 
G2 flux vacua, independent of the actual form of the Kahler potential. For nonsuper- 
symmetric vacua the results were less detailed, mainly because the constraint V — is 
quadratic in F, and, unlike the situation in [201, no regime exists in which the equations 
can be linearized. To make further progress, we will now study a class of model Kahler 
potentials for which all vacua can be computed explicitly. 

In general, at large volume, the Kahler potential is given by (EH: K = -31og(4vr 1 / 3 y x ), 
where Vx is the volume of X regarded as a function of the moduli S{. Unlike the case 
of a Calabi-Yau, where the volume function is always a third order homogeneous poly- 
nomial in the Kahler moduli, no strong constraints on Vx are known for G2 holonomy 
manifolds — just that the volume function is homogeneous of degree 7/3 and that minus 
its logarithm is convex, i.e. the second derivative of K, which gives the kinetic energies of 
the moduli, is positive definite. In general it is difficult to find simple candidate volume 
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functions which satisfy this positivity constraint. The most general homogeneous degree 
7/3 function is of the form 

n 

Vx = U s t k f^) (4-1) 



k=l 



with {at} such that 



n - 



3 

k=l 



and f{si) invariant under scaling. If we now suppose that we are in a region of moduli 
space where f(si) is approximately constant then we can take 



n 



Vx = Y[s a k \ (4.3) 



k=l 

and this in fact gives a positive moduli space metric. This justifies this particular choice 
of Kahler potentials. 

The above choice of Vx gives a simple geometry to the moduli space which is quite 
natural. The Kahler metric is 

n 3a n 3a 

ds 2 = j^dzidzi = Y 7^(^i + ds l) ( 4 - 4 ) 

8=1 1 i=l 1 

This is locally the metric of the product of n hyperbolic planes H 2 , which is: 

ds 2 = -^(dx 2 + dy 2 ) (4.5) 
l z x z 

where I is connected to the curvature tensor by 

Rn = ~l 2 ei A e 2 (4.6) 

So locally the moduli space is H 2n . Globally it is given by H 2n /Z™, because the axions ti 

are periodic variables. In this class of n-parameter Kahler potentials labeled by a«, all the 

information about X is contained in the values of the Oj. Since the moduli space metric 

(equivalently moduli kinetic terms) is singular if aj = and the moduli have the wrong 

sign kinetic terms if a« < 0, we take a» > 0. We will not restrict to any other particular 

values for the aj if it is not necessary to do so 16 . 

16 To find examples which realise these Kahler potentials, consider the case n — 7 and <Zj = 1/3, Then 
this Kahler potential correctly describes the seven radial moduli of X — T 7 and certain orbifolds thereof 
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The potential on the moduli space is given by ()3.8|) : 

c 2 / n \ 1 

V = - — ^ 3 +ya i v i sAv i s i - 3) ) + - — -=.(JV • f+ ci) 2 . (4.7) 

48vrV^V ^ / 48ttV| v ; v 7 

where v« = 

4.1 Description of the Vacua 

We now describe the vacua ie the critical points of V. The equations for the axions give 

N ■ t + c x = (4.8) 

which fixes this particular linear combination of axions. We are not concerned with 
fixing the remaining axions, since they are compact fields and will be fixed by any non- 
perturbative corrections. Our interest is in the moduli Sj. The equations of motion for the 
Si reduces to a system of n quadratic equations. For the case at hand these are equivalent 
to: 

n 

Zcijhjihj - 3) - 2/i 2 + 3hi + 9 = 0, (4.9) 

i=i 

where we defined hj = UjSj (no sum). Note that this system separates in n quadratic 
equations in one variable h{. The solutions are therefore of the form 

3 

h i = - + mi H (4.10) 

where rrii = +1 or —1, and H is determined by substituting this in (|4.9j) . This results in 
a single quadratic equation: 

5H2 -l AH - 2 ^ = °- ( 4 ' n ) 

where 

A = a-m. (4.12) 
A priori therefore, H can take two possible values: 

However, because 

HU = -H^ A) (4.14) 
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Figure 2: Left: E as a function of A = a-m. Here £(-7/3) = 3 and £(7/3) = -9/25. ffig/zi: 
dependence of A on A. At A = —7/3, A ~ 10~ 3 , and at A = 7/3, A ~ —13. The divergence at 
A = —1/3 is due to the vanishing of the volume there. 

we only ever need to consider, say, the negative branch of the square root to get all 
solutions in ([4.10)1 . In total, therefore, the number of vacua for a fixed choice of fluxes is 
2 n . We choose the following parametrisation: take all 2 n choices for m. Then 

hf> = ^+miH {ih) (4.15) 

with H = H . 

We can thus think of the vacua as the states of a system with n "spins" m,. When all 
spins are aligned with the "external field" a, that is if all rrii = +1, we have 

A = 7/3, tn=l (4.16) 
5 

When all spins are anti- aligned (rrii = — 1), this becomes 

A = -7/3, hi = 3. (4.17) 

The first of these can be shown to be the supersymmetric AdS vacuum discussed in jS] 
whilst the second is the unstable de Sitter vacuum also discussed there. The remaining 
2™ — 2 are all non-supersymmetric and could be either de Sitter or anti de Sitter. The 
metastability of these vacua will be analyzed in section I4~2l 

Substituting ()4.15|) in (|4.7|) . we get that the energy of these vacua is given by V = 
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Figure 3: Distribution of A values for = 7n/3, n = 50. 

2 

4^rE with 17 E = \E 2 - §• At fixed volume, the vacuum energy varies only through E. 
The dependence of E on A is shown on the left in fig. El However, the volume depends 
on m as well: 

^=n«r=n(pf) ni".r-. («<» 

i=l i=l Vl J|/ i=l 

so the total vacuum energy is, up to m-independent factors: 

A ~ A = E Y[ \hi\- 3a * = E |- + |- - H\~ 1+ ~. (4.19) 

i=l 

The dependence of this on A is shown on the right in fig. |21 The divergence at A = —1/3 
is due to the vanishing of the volume there, as all hi with m« = +1 vanish at this point. 
Obviously the supergravity approximation breaks down in this regime. One notable fact 
is further that the smallest positive cosmological constant is obtained when all spins are 
down, while the smallest negative cosmological constant (in absolute value) is obtained 
when all spins are up, i.e. at the supersymmetric critical point. 

At large n, the vast majority of vacua will be "halfway" the extrema. More precisely, 
if say all a« = a = 7/3n, the variable A will be binomially distributed around A = 0, as 
illustrated in El At large n this distribution asymptotes to the continuous normal density 




(4.20) 



17 The normalization is chosen such that E equals the term inside the big brackets in l|4.7|l . Also, 
4E=\{\Ff-Z\W?)- 
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For large n this is sharply peaked around A = 0, for which E = —3/20 and |F| 2 = ^c 2 ., 
so the majority of vacua are AdS and break supersymmetry at a scale M s 2 usy = |F|/V^ 2 ~ 
c 2 M^/m p . 

Thus far we have considered solutions in terms of hi. The actual values of the moduli 
are given by 

Si = -^p (4.21) 

Since the moduli fields Sj are positive in the supergravity approximation, the signs of the 
hi and flux quanta N{ must be correlated. Without loss of generality we take C2 to be 
positive. Then Ni and hi must have opposite signs. Now as long as A > —1/3, every hi is 
automatically positive, so any such vacuum must have negative Ni. When A < —1/3, hi is 
positive if m 8 = — 1 and negative if = +1, so these vacua must have sign Ni = sign mj. 
Recall that the condition A < —1/3 is also the condition to have E > 0. 

Thus, for any given rh, there is a unique choice of sign for each N{ which renders all 
s- 7 positive for all choices of |A^|. 

On the other hand, not all given, fixed fluxes N gives rise to the same number of 
vacua. The following cases can be distinguished: 

• All Ni < 0: set of vacua = {m|v4 = a ■ rh > —1/3} U {(—1, —1, • ■ • , —1)}- All vacua 
in the first set are AdS. The additional one is dS (but is unstable, as we will discuss 
in the next section). There are of order 2 n such vacua (between 2 n_1 and 2 n for 
example when all dj are equal as for the distribution of fig. EJ) . 

• Some Ni > 0, and a • sign(iV) < —1/3: just one vacuum, given by m; = sign(A r j). 
This vacuum is dS (but again will turn out to be unstable). 

• Some Ni > 0, and a • sign(iV) > —1/3: no vacua. 

As noted before, the choices of rh for which A is at or near —1/3 do not correspond 
to vacua within the region of validity of our computations for reasonable values of C2, 
because some of the moduli, and hence the volume, will be at or near zero then. 

4.2 Metastability 

We now turn to the question of metastability. We found that all dS vacua (i.e. the vacua 
with A < —1/3) for our model ensembles have a tachyon and hence are perturbatively 
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unstable. The proof can be found in appendix |B| We therefore focus on AdS vacua in 
what follows. 

In general the condition for the perturbative stability of an AdS vacuum is the 
Breitenlohner-Freedman bound [60j, which is given by 

didjV - ^V5ij > (4.22) 

where the derivatives are with respect to the canonically normalised scalars, or equiva- 
lently the Hessian is expressed in an orthonormal frame. 

The details of the stability analysis are somewhat technical and are given in appendix 
iBl We find that exponentially large numbers of the 2 n vacua are in fact metastable. 
Specifically, vacua for which A > | are always metastable. Vacua with — | < A < | 
can in principle also be metastable (and actually turn out to be local minima), but this 
is rather exceptional: they correspond to having only one of the m.j equal to +1. In 
particular, since ^ cti = 7/3, this means that the corresponding a, must be greater than 
1, and thus there cannot be more than two such solutions. 

Some lower bounds on the numbers of metastable vacua with A > 1/3 can be derived 
as follows. For simplicity, but without loss of generality, we put 

a n > a n -i > ■ ■ ■ > a\- (4.23) 

When a n > |, all solutions with m n = +1 correspond to A > ^, and we have 

iV stab > T- 1 (4.24) 

one of which is the supersymmetric solution. When a n + a„_i > | the number of stable 
vacua is at least 

iVstab > 2 n ' 2 (4.25) 

and so on. So in a model with a„+a n _ 1 + . . .+a n _ J+1 > |, but with a n + a„_i + . . . + a„_j < 
|, the number of stable vacua is at least 

iVstab > 2"- J (4.26) 

Because of f!4.23|) and the fact that di = |, we cannot have j > An/1 . So for a model 
with n moduli, the number of stable vacua is surely bigger than 

iVstab > 2 n " 4n/7 = 2 3 ™/ 7 , (4.27) 
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which is exponentially smaller than 2 n but still exponential in n. 

Actually this number is very hard to reach and for generic models the number of stable 
vacua is much bigger than this. Take for example the case aj = 3/7n of the figure, for 
which A = j A is distributed according to (J4.20)) in the large n limit. The number of 
vacua with A > 1/3 is then given by integrating the distribution (|4.20|) from A = 1/3 to 
A = oo. At large n this gives asymptotically 

iV stab /2" « exp(-n/98). (4.28) 

Again for large n this is an exponentially small fraction, but still exponentially large in 
absolute number. In fact for n = 100 the stable fraction is still about 10%. For n = 1000, 
this goes down to about 10~ 6 , but this is not a small number compared to the total 
number of vacua, which is 2 1000 ~ 10 300 . 



4.3 Distributions over moduli space 

Let us fix C2 and fh. We want to study the distributions of physical quantities over the 
space of vacua parametrized by N. As discussed in section 14.11 the sign of each is 
completely determined by m. We can therefore restrict to counting positive iV$ = |A^|. 

Let us start by finding the number of such vacua in a region TZ given by Sj > s*. By 
the substitution Sj = C2a i \h\ m * > \/Ni, this condition becomes 

ft < (4.29) 

So, in the large N approximation, the number of vacua at fixed m in this region is 

n {a) { Si > s*) = c« n ^1 = Mr n i^i vo1 ^)- ^ 4 - 3 °) 

i= i ^ \ 6 s i=i 

As in (13.22)1 . TZ is the region of the complexified moduli space projecting to 1Z. In partic- 
ular, the number of vacua in any finite region of moduli space is finite. Note that in the 
super symmetric case hi = 3/5, this reproduces the general formula 1)3.22)1 . More gener- 
ally, we see that also nonsupersymmetric vacua are distributed uniformly with respect to 
the volume form in the supergravity approximation, but that their density relative to the 
super symmetric vacua, given by EL (5^/3), is higher. Moreover, the density grows with 
increasing numbers of anti-aligned spins. The highest density is that of the dS vacua with 
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all vrii = — 1, which is 5" higher than the density of supersymmetric vacua. Obviously, 
this does not mean that in total there are 5 n times more dS maxima as supersymmetric 
critical points, since we know there is a one-to-one correspondence between them (in the 
supergravity approximation). The density in a given region is higher simply because the 
nonsupersymmetric vacua sit at larger radii. Integrated over the entire moduli space in 
the supergravity approximation, we do not run into a paradox, because both numbers 
are then infinite. In the fully quantum corrected problem, the numbers presumably will 
be finite (by analogy of what happens for type II flux vacua due to worldsheet instanton 
corrections), but then of course also the relative densities will change. 

In order to have a meaningful four dimensional effective theory, decoupled from the 
KK modes, we need the Kaluza-Klein radius to be smaller than the AdS radius. Taking 
all Si ~ s, we have 

Adb A clm? p c\m 2 p v ' 

s 2/3 3 

^kk ~ . ( 4 - 32 ) 

so R KK < Ra&s iff s > c 2 /2 . 18 

This also ensures the validity of the supergravity approximation. The number of vacua 
at fixed m satisfying this condition is given by ()4.30|) : 

n 

N(m){si > c 2 /2 ) = cl /2 JJof|/ii| (4.33) 

1=1 



Finally, to get the total density for all possible m as well, we must sum f!4.30|) over 
all m. Because of the absolute values and the nonlinear dependence of hi on m, this is 
not easy to compute analytically even in special cases. But one can get numerical results 
without much effort. For example in the case with all cij = 7/3n, we get numerically that 

n>i ~ ( 3 - 328 )" ( 4 - 34 ) 

fa i 

18 The assumption that all s, ~ s can be relaxed. Then one can prove that Rkk < Ra<is is guaranteed 

4/7 1/2 

if Si > c 2 . However for most vacua, Si > c 2 will be sufhcient to have the required scale hierarchy, so 
we stick to this estimate. 
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Figure 4: Dependence of V^ ax on A. At A = 7/3, V^ ax = (3/5) 7 / 3 « 0.3, at A = -7/3, 
ymax = 3 7/3 ra 13j and the zero i s at ^ = _i/3. 

This approximation becomes better for large n, but is quite good for smaller values as 
well. For n = 1, the exact result is 3 + 3/5 = 3.6, which is already not too far from this 
expression. 

4.4 Distributions of volumes and cosmological constants 
4.4.1 Fixed rh 

Now, let us consider the distributions of the volumes Vx- For a given in the volume is 
given by (|4.18j) : 

Vx = c 2 ^fl( a -M) aj (4.35) 

We see that as the Nj go to infinity, Vx tends to zero, so the density of vacua (strongly) 
increases with decreasing volume, i.e. large volumes are suppressed. This agrees with 
what we found earlier for the general supersymmetric case in section 13.31 
The maximal value that Vx can assume is obtained when all Nj = 1: 

\/- ax = c 7 2 /3 Yl a7\hj\ aj = c 2 /S Yl a a / V^. (4.36) 

3 3 

Here we isolated the m-dependent product Ylj \ hj\ aj : 
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The variation of Vj|? ax over different choices of rh is entirely given by the dependence of 
this function on A = a ■ rh. This is shown in fig. El 

When we take all Oj = 7/3n, and we consider say the susy case rrii = +1 so hi = 3/5, 
(I4.36j) becomes 

rj N 7/3 



VT*U=(j^J , (4-38) 

which is similar to the toy model result (J3.33j) . 

Using the usual continuum approximation for the fluxes, it is possible to get explicit 
expressions for the volume distribution of vacua at fixed rh. Taking as example again the 
case a, = 7/3n, we show in appendix IC1 that for v = Vx/V™* < 1, at fixed rh (or fixed 
A), the vacuum number density is 



dM A [v] = } 7) rr (- \nvy- 1 v-^- 1 dv. (4.39) 
(n — 1)! 

Therefore the total number of such vacua with Vx > K goes as V* 3 ™^ 7 , in agreement 
with the general estimates of section 13.31 Note that in addition here, the density of vacua 
vanishes to order {n — 1) near the cutoff v = 1. 

Thus, in accordance with general expectations, we see that large volumes are strongly 
suppressed, and more so when there are more moduli. 

Analogous considerations can be made about the distribution of cosmological constants 
at fixed rh: 

r 2 E 

A = (4 4Q) 

with E(A) as defined in section l4~Tl Clearly, the distribution of A is completely determined 
by the distribution of Vx- Thus, because large volumes are suppressed, we see that small 
cosmological constants are suppressed. In particular there is a lower bound on |A|: 

|AL in = — c 2 " 5 TT a- 3a > -J^— = — c 2 " 5 FT af a > A (4.41) 
1 1 48vr 2 A. 1 J (V^ ax ) 3 48vr 2 11 3 y ' 

with A as given by ()4.19|) and plotted in figure El 
4.4.2 All rh 

So far in this section, we have studied distributions at fixed rh, or equivalently at fixed 
A. To get complete statistics of all vacua, we need to combine these results with the 
distribution of solutions over values of A. 
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One general observation one can make is that since there are no metastable vacua at 
A < —1/3, the largest possible volume of a metastable vacuum is obtained at A = 7/3 
(and Ni — 1), the supersymmetric solution. This can be seen from fig. HJ Therefore, the 
maximal volume for a metastable vacuum is c 7 / 3 ]X/ a j J (3/5) 7 ^ 3 . It is not hard to show 
further that (7/3n) 7 / 3 < J~J • < (7/3) 7 / 3 . The former corresponds to all cij equal, the 
latter to the limiting case a, — > for all but one eij. Similar considerations hold for the 
cosmological constant. Thus we arrive at the result that for our ensembles: 

V x < 2.2 c 7 / 3 n~ 7/3 , | A| > 2.3 x KT 4 5 n 7 (4.42) 

where 1 < n < n. 

To get more refined results on the actual distributions, we need the precise distribution 
of solutions over A. This depends on the values of a«. When all a; = 7/3n for example, 
the distribution is binomial, peaked around A = 0, which for sufficiently large n is well 
approximated by the normal distribution given in (J4.2U)) . When all <2j are approximately 
equal, the distribution will still be approximately binomial, and (J4.2(J|) is still a good 
approximation for large n. When say a\ = 1 and all other a; = 4/3n with n large, the 
distribution will have two peaks, one at A = 1, corresponding to ni\ = +1, and the 
other one at A = — 1, corresponding to mi = — 1. In the following we will work with the 
distribution given in (|4.2(Jj) . 

The joint distribution for A and Vx is then given by multiplying the distributions ()4.2U|) 
and ()4.39|) . Note that ()4.39|) depends on A through V^". One can also change variables 
from A to E and thus write down a joint distribution for E and Vx, or equivalently (and 
physically more relevantly) for A and Vx, as we did for the Freund- Rubin ensemble. The 
explicit expressions are not very illuminating, so we will not get into details here. An 
example is plotted in fig. 

Clearly, Vx and A are correlated, as in the Freund-Rubin case. This follows directly 
of course from the relation A ~ E/V x . In particular given one variable, we get a roughly 
Gaussian distribution of the other variable. Qualitatively this is somewhat similar to our 
Freund-Rubin ensemble (compare 19 with fig. although the details are different. For 
example in the Freund-Rubin case, at fixed A, the vacua accumulate near the lower bound 
on Vx, whereas in the G 2 case they accumulate near the upper bound. Without constraint 
on A this gets reversed for both. 

19 Notice that fig^shows |A| _1 instead of A on the vertical axis. 
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Figure 5: Density plot of the joint distribution of AdS vacua over Vx and A, for C2 = 100, 
n = 20. 
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Figure 6: Distribution of cosmological constants dM = p(A) dA. Here we took cij = 7/3n, 
n = 50, C2 = 100, and we restricted to stable vacua with Vx > 5. 
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Finally, one can get the distribution of cosmological constants for vacua with volume 
above some cutoff value by integrating the joint density over Vx- Again for the case with 
all dj = 7/3n, we obtain the distribution for A shown in fig. El for n = 50, C2 = 100 and 
Vx > 5. The cutoff at smaller values of |A| appears because of the lower cutoff we impose 
on the volume; the lower we take the volume cutoff, the lower the value of A at which 
the density peaks. This is because 1/V X sets the scale for A. The cutoff for small |A| 
appears because the ensmemble does not contain vacua with arbitrarily large volume, as 
discussed at length before. 

4.5 Supersymmetry breaking scales 

Let us consider the distribution of the supersymmetry breaking scale, which for a fixed 
m is given by 



MLy(rn) = (y/eXgVDWiDWj 

m 

,1/2 



r ' 2 ' 7; (4.43) 



(487r) 1 /2^3/2 



where 



G , = I 1 + ?. A 2) H i + X 1 AH + 2 1 (4.44) 
m 1 3 4 I 8 64 y ' 



is a decreasing function of A, which is positive for — | < A < | and zero in the 
super symmetric solution (A = |). 

We see that the expression for M s 4 usy is the same as that of the cosmological constant, 
but with Gfh instead of E^. So the distribution of the supersymmetry breaking scale at 
fixed m is very similar to that of the cosmological constant: 

• it is completely determined by the volume distribution; 

• low supersymmetry breaking scales are suppressed (in line with the general expec- 
tations of section I3~T)1 20 ): 

• M s 2 usy has a lower bound given by 



20 although the situation here is a bit different then the case on which we focused there, namely A ~ 0, 
which cannot be obtained in the present ensemble. 
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Figure 7: Distribution of super symmetry breaking scales dM = p(M 4 )dM 4 , with M = 
M SUS y/m p . Here we took cij = 7/3n, n = 50, C2 = 100, and we restricted to stable vacua 
with Vx > 5. 

with V^ ax defined in (tOTjl . 

As for the Vx and A distributions, in order to get the complete statistic of all vacua, we 
need to consider the distribution of solutions over A, which depends on the Oj's. Similar 
to what we did for A in the previous subsection, we can compute the distribution of A for 
vacua with volume bounded by some lower cutoff. This is illustrated in fig. [7| The cutoff 
at large values of M susy appears because of the lower cutoff we impose on the volume, again 
because 1/V X sets the scale for M s 4 usy . The lower cutoff is there because large volumes are 
absent. 

5 Conclusions and Discussion 

The potential importance of the emerging ideas surrounding the landscape, e.g. for the 
notion of naturalness, is clear. These ideas should therefore be tested both experimentally 
by verifying specific predictions and theoretically within the framework of string theory. 
To make progress in the latter, one should further scrutinize the proposed ensembles of 
string theory vacua, to establish whether or not these vacua truly persist after taking 
into account the full set of subtle consistency requirements, quantum corrections, and 
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cosmological constraints [3T|- Parallel to that one should develop techniques to analyze 
large classes of (potential) vacua without actually having to go through their detailed 
constructions. Only with both efforts combined can one hope to get to a satisfactory 
analysis of the implications of the landscape picture within string theory. If nothing else, 
the latter program will teach us a lot about what kind of physics is possible in string 
theory and assist in guiding explicit model builders, as our example in the introduction 
demonstrates. 

The focus of our paper has been on the second program. To this end, we analyzed the 
statistics of Freund- Rubin and G2 flux compactifications of M theory, and compared this 
to known IIB results. 

The statistics of Freund-Rubin vacua is very different from that of flux compacti- 
fications on special holonomy manifolds. Most notably, they can have arbitrarily high 
compactification volume, and accumulate near zero cosmological constant. On the other 
hand, these vacua typically do not have a large gap between the KK scale and the four 
dimensional AdS curvature scale, so they are not really compactifications in the usual 
sense. This is no longer true j5U] if (in units in which the scalar curvature is one) the 
Einstein metric has all length scales much smaller than 1, but it remains a challenge to 
construct such manifolds with more than a modest scale hierarchy. Alternatively, one can 
imagine adding positive energy sources to lift the cosmological constant, perhaps even 
to small positive values, but no viable controlled scenario that would accomplish this is 
presently known. The problem is clear: in order to lift the cosmological constant to a 
positive value, one would have to add an energy source at or near the Kaluza-Klein scale, 
which makes 4d effective field theory unreliable and could easily destabilize the compact- 
ification. Adding such effects could also drastically change the distributions of vacua over 
parameter space. Therefore, the results we obtained for Freund-Rubin vacua should not 
be interpreted as showing there are parts of the landscape compatible with rough obser- 
vational requirements that for example strongly favor large volumes. But our results do 
show that this is possible in principle, and constructing Freund-Rubin-like vacua which 
overcome the above mentioned problems would therefore be all the more interesting. 

The statistics of G2 flux vacua has a number of universal features. One is that the 
distribution of vacua over moduli space is uniform with respect to the Kahler metric on 
moduli space. In essence, because the large volume region of moduli space is small when 
its dimension 63 is large, this implies that large compactification volumes are suppressed, 
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and strongly so if 63 is large. In particular there is an upper bound to the volume in a 
given ensemble, set by the Chern-Simons invariant c^- This is true as well for IIB flux 
ensembles studied thus far and for their presumed IIA mirror counterparts, where the 
maximal size is set by the D3 tadpole cutoff L. Since both the IIB and the G2 ensembles 
contain flux degrees of freedom that are not dual to flux degrees of freedom in the other 
ensemble, but are possibly dual to discrete geometrical deformations away from special 
holonomy, this suggests this behavior will persist when extending the ensembles to also 
sample non-flux discrete compactification data away from special holonomy. Our results 
on Freund- Rubin vacua on the other hand show that departure from special holonomy can 
produce the opposite behavior. So before general conclusions can be drawn, this question 
needs to be investigated in more general ensembles. Since the full lOd geometries of such 
extended ensembles tend to be hard to describe, an effective four dimensional field theory 
approach, as worked out in and used for example in jlTj to find M = 2 IIA flux vacua, 
would probably be the best way to make progress in this direction. 

Essentially because of the limited discrete tunability of G2 flux vacua, the scale of the 
cosmological constant is set by the volume, A ~ m^/V 3 . Since large volumes are strongly 
suppressed at large 63, it follows that small cosmological constants are strongly suppressed 
as well. This contrasts with type IIB ensembles, where the distribution of cosmological 
constants is uniform near zero. Similarly, the scale of supersymmetry breaking is set by 
rrip/V 3 ^ 2 , so small supersymmetry breaking scales are suppressed. For 6 3 not too small, 
this remains true even when taking into account tuning of Higgs mass and cosmological 
constants. F-breaking IIB flux vacua similarly favor high susy breaking scales, although 
M s 2 usy can be tuned to be small there, and the suppression was found to be independent 
of the number of moduli in the regime where M s 2 usy is much smaller than the fundamental 
scale [2111 . It is plausible that extending the G2 ensembles as discussed above would allow 
the supersymmetry breaking scale to be tuned small as well, reproducing the statistics 
of the generic ensembles of m this regime, but this would still favor higher scales. 
Another generalization which we did not consider in this work is D-term supersymmetry 
breaking. We refer to [20] for further discussion of this possibility. 

The next step to take would be to embed standard-like models in these (and other) 
flux ensembles, and study the distributions of various couplings. In particular it would 
be interesting to analyze the impact of environmental constraints on these distributions, 
to see if for example the split supersymmetry scenario of j3H] naturally emerges, and 
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if these ensembles are "friendly neighborhoods" of the landscape, in the terminology of 
[2H] • To get an idea of the possibly very important role of environmental constraints in 
determining the physically relevant distributions, consider the "atomic principle" |38j . 
which, when the light Yukawa couplings are fixed, constrains the weak scale v to lie in 
a narrow window set by the (naturally small) QCD scale. More generally, leaving the 
Yukawa couplings y arbitrary, it is actually the combination yv which is constrained to 
lie in a narrow window, so the environmental explanation of the smallness of v would 
be lost if y scans densely around zero in the ensemble. Now, in M theory vacua where 
the quarks and leptons are localized at isolated points in the extra dimensions [22 EH] , 
the Yukawa couplings are generated by M2-brane instantons wrapping 3-cycles containing 
these points and are therefore roughly of order e~ s , with s one of the geometric moduli 
we have worked with throughout our analysis. The values of s vary over the ensemble, so 
the Yukawas scan many orders of magnitude but not uniformly around zero: as we 
have seen, large values of s are suppressed and bounded above. Hence these ensembles 
are precisely of the kind envisioned in [2H] that would be sufficiently friendly to allow 
the atomic principle to play a role in explaining the smallness of v. In fact, it may even 
play a role in explaining the smallness of the light Yukawas. A priori, these ensembles 
favor smaller s and therefore larger Yukawas. However, with smaller Yukawas, the Higgs 
is allowed to vary over a wider range. Taking this into account gives an enhancement 
of order e s in the number of vacua compatible with the atomic principle, and this could 
easily overpower the a priori l/s k suppression of large s. Apart from illustrating the 
potential power of environmental principles, this also shows that such constraints can 
strongly influence the physically relevant distributions. We plan to return to these and 
other issues in future work. 
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A Normalizations 



We defined 

K = -3 In (Att^ 3 V x ) ; W flux {z) = \n,z 1 (A.l) 

K4 

The expression for the supergravity potential is the standard one for the case of dimen- 
sionless 21 scalars: 

V = k\ e K (g^DiWDjW - 3\W | 2 ) (A.2) 

Note that when expressed in terms of the moduli z, none of these effective four dimen- 
sional field theory quantities explicitly contains the fundamental scale £m, only the four 
dimensional Planck scale K4, which is the directly measurable scale in four dimensions. 
To check our normalizations, it is sufficient to verify the tension of a domain wall corre- 
sponding to an M5 brane wrapped around a supersymmetric 3-cycle S, which is Poincare 
dual to the jump in flux AG across the wall. Let's take C = for simplicity. In the 
supergravity theory, we have 



T DW = 2e K/2 \AW\ 
2 



\ANiZ % 



= T 5 vol(S), 

which is indeed the correct expression in M theory. 



B Stability analysis 

The sign of the potential for a particular solution depends on the sign of E, with E as 

defined in section |4~T1 It is negative when A = a ■ m > —1/3 and it is positive when 

A < —1/3. We discuss these cases separately. 

21 The scalars can of course be given standard dimensions by rescaling <j> = z/ka, which absorbs the k\ 
factor for the first term in (| A. 2|) . 
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B.l AdS vacua 

An AdS critical point s° does not have to be a local minimum to be perturbatively stable. 
It suffices that the eigenvalues of the Hessian of V are not too negative compared to the 
cosmological constant, and more precisely that Breitenlo liner- Freedman bound is satisfied 

d^Vis ) - ^V(s°)<% > 0, (B.l) 

i.e. this matrix should be positive definite. 

The derivatives are done with respect to the canonically normalized scalars. The 
relevant kinetic term expanded around the critical point s° is (in 4d Planck units): 

g fj d^d^ = J2^%(d„t i dn i + d^s i ) + --- (B.2) 



2 

i 

Hence the canonically normalized scalars are 



4$ 

h d ^ti + d^Si) + ■■■ (B.3) 



3ctj Si 

With this redefinition, the condition (IB.ll) becomes 



(B.4) 



2 S ° iS ° j -did 3 V{s°) - -V(s°)5 i3 > (B.5) 



3^ 



a,- a 



t 3 



where 



d idj V = ^(^(3^-2(^ + ^) + 3(^ + ^)) 

+ % (3E + 2v 2 s 2 ) S v ) (B.6) 

S i J 



r 2 

V = -^E 
48ttV| 



E = Z + ^ajVjSjivjSj -3) 



By substituting these expressions in (|B.5|) . using the variables hi = ViSi and factoring out 

c 2 

the common positive term 48 J v -i , the stability condition becomes: 

My = Ja&jiJoE - 4(/i 2 + + 6(^ + hj)) + + hi 2 ) > (B.7) 
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where 

2 
3 

Let us evaluate (jB.7J) on our solutions. 



E = -H 2 -- = — — (5 - 9A 2 + 3AV9A 2 + 15] . (B.8) 
3 8 100 V J K J 



super symmetric solution: 

^4 3 

M = — q+ — 1 (B.9) 
25 H 10 K J 

where = yfa^aj is a rank= 1 positive definite matrix. We see that M > and so 
that the supersymmetric solution is always stable. 

Other AdS (A > -1/3) solutions: 

M = Qq+Sl + Tt (B.10) 



E 4 / 3 X 2 



where 

Q = -6E, S =J + i( H + l) ' T = -AH (B.ll) 

(1 - mj) 

*<i = g ij ' ( ) 

The following relations are valid: Q > 0, T > and S + T > 0; so Qq and Tt are 
always positive definite. S has no definite sign. We have two cases: 

1. when A > |, S > 0; then 5 1 and consequently M is positive, and the corre- 
sponding solutions are perturbatively stable; 

2. when — |<A<|,S'<0;in this case the M matrix is not positive definite 
when fh has more than one entry equal to +1 and the corresponding solutions 
are not stable. Actually, if = +1 and = +1, M has the 2x2 minor 







)* : 









Q \ ah ^ ) + | I : ) (B.13) 
which has S < as one of the two egenvalues. 

On the other hand, when m has exactly one entry equal to +1, the solution is 
a local minimum and so it is stable. 

Let us prove that when A > — | and exactly one m 8 is equal to +1, the corresponding 
solution is a local minimum. 
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Without loss of generality we choose % — 1. crV is proportional to M = q 



Q Q 



t, 



where S" = S + f-E is negative for all considered values of A. Let us call B = > 0, 



C = -%>0,Vi 



RTi. So \v | 2 = 7/3, f i > 1 and 



M = vv t + 



V 



B 



\ 



B ) 



;b.i4) 



By direct study of its eigenvalues, it can be shown that this matrix is positive definite for 
n = 2. Now we prove by induction that this matrix is positive definite for each n. 
Let us assume that G W l 



tfyiny = (v m +v-yY + By z - Cyf > 



(B.15) 



where V = (Vi,V). Let u be the extension of v E W l to 
\/x e R n+1 there exists y e W a , such that 



, i.e. ui 



Then 



Vi—x\ y ■ v = x ■ u y 2 = x 2 



(B.16) 



It follows that 



x M n+1 x = (u 1 x 1 + u • x) + Bx 2 - Cxf > 



(B.17) 



QED. 



B.2 dS vacua 

The solutions with positive potential (i.e. A < —1/3) must be local minima to be 
metastable, and so the matrix M' = Qq + S' 1 + Tt (in the notations introduced in 
the first part of this appendix) must be positive definite. We now show that this is never 
the case. 

We have q = vv t , where V{ = y/al and so |t>| 2 = 7/3. Let us consider the n x n matrix 
defined by M" = Q q + S' 1 + T 1. If M' is positive definite then M" is positive too. 
Clearly, its eigenvalues are S' + T, which is positive for A < —1/3, and 

\ = Q\v\ 2 + S' + T= 7 -Q + S' + T. (B.18) 
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By studying A as a function of A, one finds that it is positive for A > —17/21. So M" is 
not positive definite for A < —17/21 « —0.8095 and so M' is not positive for the same 
values. 

Moreover one can study the k x k submatrix of M' obtained by restriction to the 
subspace on which rrij = +1 (k is the number of such m/s). This is equal to Qvv t + S' lk, 
where we are considering only Vj corresponding to rrij = +1. Its eigenvalues are S', which 
is positive for A < —1/3, and 

^ = QJ2 + S' = Q(l + ^) + S>. (B.19) 

i=l 

By studing \x as a function of A, one finds that it is negative for a < A < — ~, where 
a ~ —1.44075. For these values of A we have a nonpositive minor and so M' is not 
positive definite. 

By combining these two results, we see that for no value of A < — | the corresponding 
solution is a local minimum. Thus, all dS critical points are unstable in our model 
ensembles. 



C Volume distribution 



Here we give details on how to arrive at the distribution (|4.39|) of volumes v = Vx/V : 
when all cii = a = 7/3n: 

dM[v]/dv = "£6(v-l[Nr) 

Ni>l » 

w / d n NS(v -T[N- a 

JNi>l 



max 
X 



.1 

\nv \n— 1 



- v- {1+1/a) \ a \, 0(1 - v) 
a [n — I)' 

(3n\ n 



n-1 ! 
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